Many-body effects in iron pnictides and chalcogenides — 
non-local vs dynamic origin of effective masses 
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We apply the quasi-particle self-consistent GH^(QSGW) approximation to some of the iron pnic- 
tide and chalcogenide superconductors. We compute Fermi surfaces and density of states, and 
find excellent agreement with experiment, substantially improving over standard band-structure 
methods. Analyzing the QSGVK self-energy we discuss non-local and dynamic contributions to ef- 
fective masses. We present evidence that the two contributions are mostly separable, since the 
quasi-particle weight is found to be essentially independent of momentum. The main effect of non 
locality is captured by the static but non-local QSG Ineffective potential. Moreover, these non-local 
self-energy corrections, absent in e.g. dynamical mean field theory (DMFT), can be relatively large. 
We show, on the other hand, that QSGVFonly partially accounts for dynamic renormalizations at 
low energies. These findings suggest that QSGVF combined with DMFT will capture most of the 
many-body physics in the iron pnictides and chalcogenides. 

PACS numbers: 74.70.Xa,71.27.-|-a,71.15.-m 
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The discovery of superconductivity in the iron pnictides 
and chalcogenides has triggered much effort into under- 
standing their electronic properties, [ij The first theoreti- 
cal insight into the pnictides was gained by density func- 
tional theory (DFT) within the local density approxima- 
tion (LDA), which correctly predicted the Fermi surfaces 
of LaFeP0[2[ and LaFeAsO[3[, as well as the striped an- 
tiferromagnetic spin ground stateQ. However there is 
ample experimental evidence for sizable correlation ef- 
fects in the pnictides, manifesting themselves in high ef- 
fective masses as witnessed by optical spectroscopy [5[, de 
Haas-van Alphcn measurements @, or bandwidth renor- 
malizations from photoemission spectrosco py (P ES) 0- 
lll, as well as low magnetic moments [ij, ll3[. These 
issues have been successfully addressed by many-body 
techniques, yielding correct effective masses Il4 |l5[, or- 
dered moments [13], and good structures flil Il7|. 

In many of these works, however, correlation effects 
have been accounted for by treating local interactions 
for a subspace of orbitals only. While this improves the 
description, several shortcomings persist : (a) correct- 
ing only a subspace means that large parts of the elec- 
tronic structure are still DFT derived. The choice of 
the exchange-correlation functional being discretionary, 
moreover causes a dependence on the effective one- 
particle starting point. Also, results depend on the 
correlated subspace chosen, (b) Out-of-subspacc sclf- 
energies[ll| are neglected, leading e.g. to underestimates 
in p-d gaps, (c) Non local interactions are treated on the 
effective one-particle level. To address these issues, we 
applied the quasi-particle self-consistent G approxima- 
e iron pnictides and chalcogenides. We 
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discuss band-structures, Fermi surfaces, and density of 
states (DOS) and elucidate the origin of many-body effec- 
tive masses. The essence of our findings is : (1) Non-local 
self-energy effects are not small. These are neglected in 



DMFT[22[ based approaches. (2) QSGW^, which includes 
non-local correlations, produces excellent Fermi surfaces, 
but (3) does not adequately account for the dynamics of 
the self-energy. (4) Non-locality and energy dependence 
are shown to be essentially separable. (5) Together, this 
indicates that a combined QSGW+DMFT approach is a 
very promising avenue in electronic structure theory. 

An important goal of the many-body theory of 
solids is to capture the one-particle Greens-function 
G, which is written as G^^ = w - iJ°(k) - S(k,a;), 
with H'^ a reference one-particle Hamiltonian, and S 
a self-energy that is defined with respect to correla- 
tion effects already included 
pie, in LDA-hDMFTjii, H^- 

„,LDA 



in 



4^ . For exam- 

-V^ -I- Vcrystal + 

VHartree + ^^Jc": ^hc Kohu-Sham Hamiltouian, with 
^={TP"^^{lo) ~ Edc)LL'RR'\RL){R'L'\. The self-energy 
j^DMFT g^j,^g ^ gg|- q£ correlated orbitals \RL), 
and the portion Edc of already contained in 

must be subtracted. In the QSGl^F approxima- 

tion, H''=H^^^^=-V^+V,rystal+VHartree+v'^l''''\ and 

Y^QSGW^QQSGW^r _ yQSGW^ ^J^gj.g ^J^g (gtatic) QSGW^ 

exchange-correlation potential, v'^l°^ , is chosen so as to 
closely mimic the same quasi-particles E\^j in TJQsgw 
in (w - TJQSGw _ SQSGH.)^ These arc given by 

[i7«^''-(k) -f 5ftS«^-^-(k, E^,)] \^^,) = i?k,]*k,) (1) 
and it can be shown [lij that a good choice for 

\ ^ [S«^''"'(k,i?k.) + S«^°^-(k,i?k,)] (*k,|. 

ijk 

(2) 

Results for paramagnetic BaFe2As2 are shown in 
Fig- HI As is clearly visible, QSGFK substantially narrows 
the iron 3d bands (by 16%) relative to the Kohn-Sham 
spectrum within LDA [4^. This narrowing corresponds 
to an enhancement of the effective mass, m'^^^''^ /m^^^, 
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FIG. 1: paramagnetic BaFe2As2. (a) bandstructures within LDA and QSGH^ (notice the Fe-3d bandwidth narrowing as 
indicated by the arrows), (b) comparison with ARPES[2^ around the F point. The slight constant energy shift between theory 
and experiment possibly originates from doping in the measured sample, BaFei.gsCoo.is As2 [13]. 



and is slightly larger than the 12% found for elemen- 
tal iron[20| (See Tab. H] for other compounds). In 
Fig. [TJb) we compare QSGW and LDA bands near the 
Fermi level Ep and k near zero, to angle resolved pho- 
toemission (ARPES) experiments [2^. It is consensus 
from experimental data that there is, at the T point, 
an outer hole pocket, while closer to T there are two 
xz/yz excitations 0, [23l . [23 . The top of these inner 
bands is very close to zerol23|, particularly for Co-doped 
samples (23. 25 1. The character of the outer pocket is de- 



bated : a mixture of xy and x^-y^ l23|, z'^ additions [26[, 
and a change of character along fc2[26[ is reported. The 
Fermi vector has been measured for the outer band only : 
/cF^O.OT-O.irA-iQlillill. QSGI^hardly effects this 
pocket; and it indeed is mainly of xy character for both 
LDA and QSGW. The Fermi wave vector is fcF=:0.24A~^ 
The inner pockets are of xz/yz character, in accord with 
experiment [i^, [2^. Relative to LDA, their size is re- 
duced to fcF=0.05A~i fc_F=0.1lA" 
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While QSGFF reduces the group velocity relative to 
LDA, the dispersion at kp near F as measured by ARPES 
is still a factor of two or so smaller. This suggests that 
QSGM^does not fully account for all many-body renor- 
malizations in this compound, as we detail below. 

We next consider the chalcogenide FeSe (Fig. [5]) . The 
QSGW band-structure displays a remarkable narrowing 
of the iron 3d bands of 22%, equally visible in the DOS. 
Also higher energy bands get renormalized : Se-4p ex- 
citations are pushed down in energy, notably improv- 
ing the agreement with PES.[llj In Fig. [^Jc-e) we com- 
pare the Fermi surfaces of FeSe and FeTe within QSGl^F 
and LDA with measurements of the Fe1.04Teo.66Seo.34 
alloy 28j . In experiment there are three hole pockets 
at F and two electron pockets around M. QSGM^and 
LDA predict similar surfaces for the xy pocket labeled 
"7", and the Fermi wave vectors (FeSe : kj, = 0.29A~-^, 



FeTe : kj, = 0.43A-1) encompass the ARPES value for 
Fe1.04Teo.66Seo.34 {kl = O.3A-I). The QSGW^ Fermi 
vectors of the xz/yz "a" and pockets shrink rela- 
tive to the LDA, and become more anisotropic. Indeed, 
the "a" pocket has kf = OMA'^ on the FA/ line for 
FeSe, and completely disappears in FeTe. The magnitude 
agrees well with ARPES measurements {kp = 0.03A~-'^), 
whereas the ARPES (3 pocket (fc^ = O.Uk"^) is some- 
what smaller than in QSGW^(fc^ = O.I9A-1 and 0.23A-1 
in FeSe and FeTe, respectively). 

For the Ill-family, we show in Fig.[3]the Fermi surface 
of LiFeAs within QSGW^, LDA and ARPES[i3. QSGW^ 
and ARPES agree very well for both large and small 
pockets: note in particular how QSGW shrinks the in- 
ner pockets at F. Both calculations and experiment [2^ 
concur in the orbital characters along TX : xy, yz, and 
xz, for the outer, middle and inner band, respectively. 
Within QSGW^, the latter has moved below Ep. At M 
the smaller electron pocket is of mixed xz/yz; the big- 
ger one, of xy character, is too small in LDA. It be- 
comes larger in QSGW, as was previously found within 
LDA-f-DMFTE3 



We turn to the mass enhancement relative to the LDA 
band masses, which have been used as a reference in the 
analysis of experiments. The enhancement is given as the 
ratio of the magnitude of the LDA and the QSGW group 
velocities near the Fermi level. The latter are 



^ (^kd^fc, (g'^=^°'"+ ^S]Q^"-"-(c^ = 0)) 

rffca (1 - (*i„|a„Ki]«^°»'|«'k.)):lo 



(3) 



A change in the velocity, and thus the effective 
mass, is possible through (a) the dynamical part of 
the self-energy via the quasi particle weight Zk = 
1/(1 -5,^3?S«^''"'(fc,w))^^o. Noting that i/Q^^"' 4- 

5^QSG«.^_y2 ^ ^^^^^^^^ ^ VHartree " GQ^'^"'^ thcrC is (b) 
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(a) band-structure of FeSe (b) FeSe density of states vs. photoemission 




(c) FeSe (d) Fc1.04Teo.66Seo.34 (e) FeTe 

FIG. 2: paramagnetic FeSe. (a) band-structure within LDA and QS GW (notice the substantial bandwidth narrowing 
indicated by the arrows), (b) DOS in comparison to photoemission [Til] . Also shown is the spectral function A{ij) of the QSGW 
that takes into account lifetime effects. Arrows indicate the Fe-3d bandwidth; grey lines are guides to the eye. (c)-(e) Fermi 
surfaces : (c)/(e) FeSe/FeTe within QSGW and LDA, (d) experimental Fermi surface of Fe1.04Teo.66Seo.34 from ARPES[2^. 
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TABLE I: effective masses and quasi-particle weights. 

^Qs&wy^LDA (jgjjQ^gg ^jjg ratio of the iron 3d bandwidth 
within QSGW^and LDA. The quasi-particle weights Z'^^^"' 
are extracted from S'S^'^* at the F point in the band basis. 
DMFT masses from Ref. [l3 . except for the second values of 
LiFeAs (Ref. [13) and FeSe (Ref.[l|). ARPES effective masses 
are obtained with respect to LDA, in accordance with Eq. 



a renormalization of the velocity through non-local corre- 
lations as encoded in G''^^°^W, and (c) an effect through 
a change in charge density. 

In DMFT approaches, the self-energy is local by con- 
struction, and enhanced masses emerge solely through 
the energy dependence of the self-energy (indeed 
m/m°"''^=2'°"^^). While in many correlated materi- 
als, in which the physics is controlled by the on-site 
Coulomb interaction, there is evidence that this is a good 
approximation[3^ . it is not clear a priori whether such an 
ansatz is warranted for the pnictidcs and chalcogenides. 



QSGW, on the other hand, accounts, albeit perturba- 
tively, for all mechanisms of effective masses : Indeed, 
the momentum-dependence is manifestly included via 
Eq. ([2|), and the energy slope enters through the QSGW 
procedure, Eq. ([!]), that determines £'k. By construe- 
tion, Eq. (O also reads ^ = {^k-Mdk^H^''"''') 
hence wjj?'''" accounts for, both, the dynamic and non- 
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In order to discuss non-local and dynamic contribu- 
tions to effective masses, we analyze the QSGW self- 
energy : Tab. [J summarizes our results for the dy- 
namical renormalization factors l/Z"^^"^', the net 3d 
bandwidth-narrowing, TO*^^'^"'/m'"°*[4^, as well as re- 
sults from DMFT and ARPES. This reveals : (a) The 
inverse quasi-particle weights are larger than the actual 
bandwidth narrowing by a factor of two or more. One 
clear source of the discrepancy originates in the nonlocal- 
ity of the self-energy, missing in both LDA and DMFT. 
Nonlocality tends to delocalize g uasi-particles [ssj . see 
Eq. ((3|) , which partially cancels the effects of the energy 
dependence of S. We analyze this further below, (b) The 
band-width narrowing m'^^"^ /m^^^ is too small with re- 
spect to experiment. This is generally expected, since the 
reduction of the quasi-particle weight Z in correlated sys- 
tems is a largely non-perturbative effect - a realm where 
DMFT excels. In particular, a DMFT studyQ ratio- 
nalized that the large effective mass in FeTe is owing to 
the local Hund's rule coupling, a multiplet effect not well 
treated in GW. Despite the insufficient dynamic renor- 
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malization, however, the trends in Z along the series of 
materials are captured, including the orbital differenti- 
ation. Indeed hole-like excitations of xy character are 
heavier than xz/yz ones for all compounds considered. 

To make a direct connection with DMFT, we 
introduce a local basis set (indexed V) by con- 
structing maximally localized Wannier functions [s^. 
35| for the iron 3d and pnictogen/chalcogen 4p or- 
bitals. We further introduce the momentum variance 

(4^)^=l/^L^EkLL' \^^\v - J^^ii'T °f a quantity 
X with respect to its local reference ^|°£/=X] 
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To quantify static corrections in QS GW beyond LDA 
at the Fermi level we define the many-body correction 
^ ^ ^QSGw _^ fjQSGw _ ijLDA^^^j^^ computc Its vari- 

ance 5fcS(a;=0) [H]. We find 4S=0.08eV for BaFeaAsa, 
and 5fcI]=0.09eV for FeSe, while the variations of the 
LDA exchange-correlation potential, are 0.22eV 

and 0.2eV, respectively. Hence, QSGW^ overall reduces 
the momentum dependence, consistent with the observed 
bandwidth-narrowing, see Eq. ([3]). The relative change 
in the momentum variation, JfeS/^fcwJ;^?*, is important in 
both compounds : 36% for BaFe2As2 and 45% in FeSe. 
This shows that the assumption of a purely local self- 
energy correction (a la DMFT) is not fully warranted. 

To quantify the momentum dependence of the quasi- 
particle dynamics, we introduce the variance A^Z = 



Z'fll'^ of the Fermi liquid weight. We find 



AfcZ«0.5% for all calculations performed here|49|. 
This establishes that non-local and dynamical correla- 
tion effects are essentially separable. Thus the self-energy 
is describable by the ansatz I](k, w) F(k) -f (1 — 
Z~^)u! within the Fermi liquid regime, which encom- 
passes (within qSGW) the range |a;| < 2eV. The QSGW 
self-consistency further imposes : F(k) = Z~^Ey^ — ej^°*. 
To introduce non-local correlations to DMFT, or. 



reversely, to improve on local fluctuations in GW, a 
GW+DMFT scheme has been proposedjsl]. In the sim- 
plest of implementations local and non-local self-energies 
from DMFT and GW, respectively, are combined. This 
was rationalized further by noting that self-energy contri- 
butions beyond GFFare mostly localjlll, thus amenable 
to DMFT. Here, we showed that even on the GW 
level, the (low-energy) dynamics is essentially local, too. 
Hence, non-local correlations can be accounted for by the 
static QSGW potential. This has several advantages : 
Due to its complexity, GFK-I-DMFT has, so far, only 
been employed non-selfconsistentlv(36i [Sq. This defies 
many of the method's virtues, e.g. it persists a start- 
ing point dependence, and issues of double counting. 
Here, we advocate to circumvent this by imposing (quasi- 
particle) self-consistency for the GW^only, i.e. introduc- 
ing a QSGFK-I-DMFT scheme. The latter is manifestly 
void of any DFT dependence, and the double counting is 
well defined [50l |. 

For the iron pnictides and chalcogenides we conclude 
that the QSGW approach yields very good Fermi sur- 
faces, as well as sizable corrections to spectral features 
at finite energies. We further showed that non-local cor- 
relations are important, and can be incorporated through 
the QSGl^F effective potential. Moreover we showed that 
dynamical correlations are mostly local. They are not 
sufficiently accounted for in QSGW, especially when cor- 
relations are strong, as in the chalcogenides. These find- 
ings suggest that QSGW and DMFT taken together will 
be sufficient to incorporate most of the correlations in 
electronic structure. 

The support of the DOE-CMSN grant "Computational 
Design of Fe Based Superconductors" is gratefully ac- 
knowledged. 




FIG. 3: Fermi surface of LiFeAs. = 0-plane in the Bril- 
louin zone for 2 Fe atoms; experimental intensity from Ref . [lol. 
Notice how QS G drastically shrinks the inner pockets at F. 
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A matrix structure in the space of orbitals is assumed. 
Our LDA results (method of Ref. Is^ ) are congruent with 
previous works, e.g. Ref. 

Interestingly, the spin-orbit coupling makes a qualitative 
difference here. Indeed without it, the pockets are re- 
versed in LDA, thus yielding wrong orbital characters. 
For insulating VO2 correlation effects from a DMFT clus- 
ter extension were found to be mappable onto a non-local 
one-particle potential^ljl akin to the QSGH^ procedure. 
We find the ratio of bandwidths is very similar to the 
inverse ratio of group velocities on the Fermi surface. 
We restrict L, L' to the 3d orbitals of iron. 
The Wannier functions based on the LDA and the QS G"W^ 
are of course different. However, as the Wannier spread 
differ by less than 3%, we neglect this effect. 
Afc [1 - (9„KE«'^'^'"'(fc,a;)]"^<Z/10 for |a;|<2eV : the mo- 
mentum variance of the quasi-particle dynamics is small, 
namely by the local projection of G'''^'^'^W. 



